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ABSTRACT 

This  paper  first  generalizes  a characterization  of  polyhedral  sets 
having  least  elements,  which  is  obtained  by  Cottle  and  Veinott  [6] , to  the 
situation  which  Euclidean  space  is  partially  ordered  by  some  general  cone 
ordering  (rather  than  the  usual  ordering) . We  then  use  this  generalization 
to  establish  the  following  characterization  of  the  class  Q of  matrices 
(C  arises  as  a generalization  of  the  class  of  Z-matrices,  see  [4],  [13], 
[14]):  M e C if  and  only  if  for  every  vector  q for  which  the  linear  com- 


plementarity problem  (q,M)  is  feasible,  the  problem  (q,M)  has  a solution 
which  is  the  least  element  of  the  feasible  set  of  (q,M)  with  respect  to  a 
cone  ordering  induced  by  some  simplicial  cone.  This  latter  result  generalizes 
the  characterizations  of  K-  and  Z-matrices  obtained  by  Cottle  and  Veinott 
[6]  and  Tamir  [21]  respectively. 

AMS(MOS)  Subject  Classification  - 90C99,  06A10,  15A39,  15A48 

Key  Words  - Cone  orderings,  least  elements,  linear  complementarity,  generalized 
Leontief . 

Wor)t  Unit  Number  5 - Mathematical  Programming  and  Operations  Research 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024  and  the 
National  Science  Foundation  under  Grant  No.  MCS75-17385. 


EXPLANATION 


The  linear  complementarity  problem  is  nowadays  an  extremely  important 
subject  in  mathematical  programming.  Its  wide  applications  can  be  found  in 
such  diverse  areas  as  optimization,  economics,  structural  engineering,  finance, 
free  boundary  value  problems  and  optimal  stopping.  The  area  of  large-scale 
linear  complementarity  problems  has  in  recent  years  received  an  increasing 
amount  of  interest  because  of  its  potential  usefulness  in  solving  discretized 
partial  differential  equations.  A fundamental  problem  arising  here  is  the 
solution  of  these  large-scale  problems.  Many  of  the  currently  available  linear 
complementarity  algorithms  are  either  not  applicable  or  not  suitable  to  solve 
large-scale  problems  because  of  limited  computer  storage  and  expenses. 

Recently,  Professor  Mangasarian  has  proposed  solving  a class  of  large-scale 
linear  complementarity  problems  as  linear  programs.  The  purpose  of  this  paper 
is  to  provide  a characterization  of  this  class  of  problems  in  terms  of  some 
geometric  notions  in  Euclidean  space. 
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ON  CONE  ORDERINGS  AND  THE  LINEAR 
COMPLEMENTARITY  PROBLEM 

Jong-Shi  Pang 


1.  INTRODUCTION 

Recently,  there  are  a number  of  papers  in  the  literature  which  are  concerned  with 
characterizing  polyhedral  sets  having  least  elements  [61 , with  characterizing  certain 
classes  of  matrices  in  terms  of  linear  complementarity  problems  having  least-element  solu- 
tions [6] , [16] , [21] , and  with  solving  linear  complementarity  problems  as  linear  programs 
[41,  [5],  [131,  [141,  [151,  [181.  In  fact,  these  three  subjects  are  very  closely  related 
to  each  other.  Among  those  papers  mentioned  above,  the  first  one  [6]  seems  to  be  the  prime 
motivation  for  investigating  the  various  relationships  botv.’cen  the  three  theories.  The 
essential  result  obtained  in  that  paper  is  a theorem  whiclj  characterizes  polyhedral  sets 
having  least  elements  with  respect  to  the  usual  ordering  of  Euclidean  space.  As  an  ap- 
plication of  this  characterization,  the  authors  of  the  paper  derived  a characterization  of 
the  class  of  K-matrices  in  terms  of  linear  complementarity  problems  having  least-element 

solutions.  It  thus  follows  that  linear  complementarity  problems  with  K-matrices  can  be 

* 

solved  as  linear  programs  . The  characterization  of  K-matrices  was  later  extended  to  Z- 
matrices  by  Tamir  [21) . Therefore,  linear  complementarity  problems  with  Z-matrices  can 
also  be  solved  as  linear  programs. 

An  application  exploiting  the  fact  that  linear  complementarity  problems  witli  Z- 
matrices  have  least-element  solutions  has  been  described  in  [16).  Sec  also  [17).  In  [16), 
the  author  studied  a class  of  large-scale  linear  complementarity  problems  arising  from 
quadratic  programs  with  upper  and  lower  bounds  on  the  variables  and  with  no  other  con- 
straints. These  quadratic  programs  have  many  aiiplications  in  various  areas.  See  [3). 

He  [16]  presented  a fast  and  efficient  algorithm  for  this  class  of  linear  complementarity 


* See  [4]  for  a more  precise  description  of  this  concept.  Througliout  the  paper,  the 

phrase  "solving  linear  complementarity  problems  as  linear  programs"  has  the  meaning  as 
described  in  [4] , I13J , 114) . 
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problems  which  makes  use  of  an  efficient  algorithm  (2]  to  find  the  least-element  solution 
of  a linear  complementarity  problem  with  a Z-matrix. 

Mangasarian,  in  two  recent  papers  [13],  [14],  introduced  certain  new  classes  of 
matrices  for  which  he  showed  that  the  corresponding  linear  complementarity  problems  can  be 
solved  as  linear  programs.  His  method  of  derivation  has  nothing  to  do  with  least  elements. 

In  [14],  R.  w.  Cottle  and  the  author  summarized  the  results  in  Mangasarian' s two  papers  by 
showing  that  all  the  classes  of  matrices  studied  in  the  paper  - that  are  seemingly  quite  dif- 
ferent - are  in  fact,  subclasses  of  a large  class  of  matrices,  which  we  have  denoted  by  C • 
(Incidentally,  C includes  the  classes  of  Z-  and  K-matrices.)  Moreover,  we  have  shown  that 
Mangasari an' s results  can  indeed  be  derived  via  Cottle  and  Veinott's  theory  of  polyhedral 
sets  having  least  elements,  thereby  tightening  up  the  connection  between  this  latter  theory 
and  that  of  solving  linear  complementarity  problems  as  linear  progriims . 

It  is  natural  to  ask  whether  matrices  in  C can  in  fact  be  characterized  in  terms  of 
linear  complementarity  problems  having  least-element  solutions.  That  tliis  might  be  pos- 
sible is  suggested  by  the  fact  that  it  is  possible  for  tlie  classes  of  Z-  and  K-matrices. 

Our  purpose  in  this  paper  is  to  provide  a positive  answer  to  this  question.  Wo  would  like 
to  mentitm  that  a generalization  of  c has  been  studied  in  [5]  (least-element  aspect)  and  in 

[15]  (non-least-elc  ment  aspect).  Moreover,  it  has  been  shown  [If!]  tliat  C is  closely  re- 
lated to  the  well-known  class  Y of  matrices  wliose  cbaract eri zati on  has  long  been  an  open 
problem  in  the  theory  of  (he  linear  complementarity  juoblcm  but  has  recently  been  established 
in  [18]. 

In  order  to  characterize  (3  , we  need  to  consider  Euclidean  space  as  being  partially 
ordered  by  some  general  partial  orderings  (rather  than  just  the  usual  ordering  that  is 
always  implicitly  implied  in  all  the  known  characterizations).  These  partial  orderings 
are  induced  by  pointed  cones  and  are  thus  called  cone  orderings.  They  certainly  include 
the  usual  ordering  as  a special  case.  We  shall  develop  a theory  of  polyhedral  sets  having 
elements  that  are  least  with  respect  to  these  cone  orderings.  The  theory  is  an  extension 
of  that  obtained  in  [6] . The  key  characterization  theorem  is  described  in  terms  of  a gen- 
eralized Lcontief  property  presented  in  .Raigal  [19]  and  will  lx:  used  to  characterize  C 
the  manner  described  in  the  last  paragraph.  As  in  (<0  , we  do  not  address  the  question  of 
the  existence  of  the  least  elements.  Instead,  we  refer  the  interested  reader  to  [17]. 


The  plan  of  this  paper  is  the  following.  The  next  section  is  a summary  of  back- 
ground materials.  It  contains  two  parts.  In  the  first  part,  we  review  some  basic  defini- 
tions and  fix  our  notations.  In  the  second  part,  we  state  some  known  results  that  are  im- 
portant to  the  development  of  our  theory.  In  the  third  section,  we  consider  r"  as  being 
partially  ordered  by  some  cone  ordering  and  develop  a theory  of  polyhedral  sets  having 
least  elements.  We  include  a theorem  characterizing  a (strictly)  isotone  linear  function 
on  r"  under  a cone  ordering.  This  theorem  is  believed  to  be  new  and  is  related  to  the 
problem  of  finding  the  least  element,  provided  that  it  exists.  In  the  fourth  and  last 
section,  we  establish  the  promised  characterization  of  (J  . 
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2.  BACKGROUND 


2.1.  Basic  Definitions  and  Notationn.  Throughout  this  paper,  will  denote  the  non- 

negative orthant  of  Euclidean  n-space  r",  and  r"^^  will  denote  the  class  of  real  nxS 

matrices.  We  denote  the  i-th  column  (row)  of  a matrix  A e r''^^  by  A . (A.  ).  By  e^ 

•IX" 

we  denote  the  i-th  unit  vector,  i.e.  the  vector  whose  components  are  all  zero  except  the 
i-th  component  which  is  one.  If  Sc  r",  we  denote  the  interior  of  S by  int(S)  . 

For  a given  vector  q c r"  and  matrix  M e r"^",  wo  denote  the  linear  complementarity 
problem  of  finding  a vector  x e r"  satisfying 


q + Mx  ^ 0 and  x (q  + Mx)  = 0 

by  the  pair  (q,M)  . Note  that  the  nonnegativity  of  the  vectors  x and  (q  ^ Mx)  are 
meant  comiionentwise.  By  the  feasible  set  for  the  problem  (q,M)  , we  mean  the  polyhedral 


X(q,M)  = {x  f R^:  q + Mx  O)  . 

The  problem  (q,M)  is  said  to  be  f easible  if  X(g,M)  ^ <()  . 

A matrix  M e.  R^^  is  said  to  be  a Z-matri  x if  it  has  nonpositive  off-diagonal 

entries.  The  Z-matrix  M c r'^'^  is  said  to  be  a K-matrix  if  there  is  a vector  x t R^ 

+ 

such  that  Mx  > 0 (componentwise).  Properties  of  these  two  classes  of  matrices  )iave  been 

nXH 

surveyed  in  (8).  Let  A c R . It  is  said  to  be  Leontief  if  it  has  at  most  one  positive 

I 

element  in  each  column  and  there  is  a vector  x e R^  such  that  Ax  > 0 (comijonentwise)  . 


If  A is  Leontief,  then  the  system 


Ax-b>0,  x>0 


is  called  a Leontief  substitution  system.  Properties  and  applications  of  Leontief 

matrices  and  of  the  associated  Leontief  substitution  systems  arc  well  recognized  in  the 

literature.  See  (7),  [22J  . It  is  clear  that  K-matricos  are  Leontief;  and  conversely, 

nX?' 

square  Leontief  matrices  with  positive  diagonal  elements  are  K-matrices.  For  A e R 
with  full  row  ran)c,  we  say  an  nX  n submatrix  B is  a basis  if  it  is  nonsingular. 

A partial  ordering  S (see  [1))  on  a set  S is  a binary  relation  on  S,  which 
satisfies  for  all  x,y  and  z c S,  the  following  three  axioms: 


(PI)  X < X (ref lexivity) ; 
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(P2)  X < y and  y < z imply  x < z (transitivity); 

(P3)  X < y and  y < x imply  x = y (antisymmetry) . 

The  set  S is  said  to  be  partially  ordered  (by  ^ ) or  a poset  if  <,  is  a partial  ordering 
on  S . We  denote  the  poset  S together  with  the  partial  ordering  < by  the  pair  (S,  ^ ) 
It  is  clear  that  every  subset  of  a poset  is  a poset  with  the  induced  ordering.  If  x and 
y are  elements  in  a poset  S,  x < y and  x ^ y,  then  we  write  x < y . 

Example . The  usual  ordering  ^ of  r"  is  defined  as  follows:  For  x,y  e R^,  x < y if 
and  only  if  x^  ^ y^  for  every  i . It  is  trivial  to  show  that  this  is  a partial  ordering 
Later  in  this  section,  the  usual  ordering  of  r"  will  t>e  generalized.  For  this  particular 
ordering,  we  write  x < y to  mean  x^  < y^  for  every  i . This  is  not  to  be  confused 

with  < which  is  used  for  other  orderings  and  has  a weaker  meaning. 

* 

Let  T be  a subset  of  the  poset  (S,  < ) . An  element  t e T is  a least  element 

* 

of  T (with  respect  to  or  under  < ) if  t ^ t for  every  t c T . The  least  element  of 
a poset,  if  it  exists,  must  bo  unique.  This  follows  immediately  from  its  definition  and 
the  antisymmetry  of  the  partial  ordering. 

Let  (S,  < ) and  (S',  <')  be  posets.  A mapping  f:  S S'  is  said  to  be  i sotono 
if  x,y  c S and  x < y imply  f(x)  S'  f (y)  . An  isotone  mapping  is  strict  if  x<  y 
implies  f(x)  <'  f(y)  . When  S'  is  the  real  line  and  S is  the  usual  ordering  ^ of 
scalars,  we  say  that  f is  a (strictly)  isotone  real -valued  function  if  the  mapping 
f:  (S,  < ) (R,  ^)  is  (strictly)  isotone. 

We  review  a few  concepts  about  cones  in  Euclidean  space  [20] . A subset  C of  r" 
is  called  a cone  if  it  satisfies  the  following  three  conditions: 

(Cl)  0 e C; 

(C2)  Xx  e C for  every  X c R^  and  x e C; 

(C3)  X + y £ C for  x,y  c C . 

The  cone  C is  said  to  be  pointed  if  x £ C and  -x  £ C imply  x = 0.  It  is  finitely 
generated  if  there  exists  G £ such  that  C = {x  e r":  x •=  Gy  for  some  y £ R™}  . 

In  this  case,  we  denote  C by  pos(G)  . A cone  C is  polyhedral  if  there  exists  F£  r’’^" 
such  that  C = {x  £ r”:  Fx  ^ o}  . A cone  C is  simplicial  if  C » pos (X) 
for  some  X £ r"^"  and  X is  nonsingular. 
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Note  that  conditions  (C2)  and  (C3)  together  imply  that  cones  are  convex.  It  is  well- 
known  (see  [20])  that  a cone  is  finitely  generated  if  and  only  if  it  is  polyhedral.  For 
an  arbitrary  subset  S c r",  we  define  its  polar  cone 

* n T 

S ={ycR:  xy^O  for  every  x c s}. 

* 

Obviously,  S is  a nonempty  closed  cone.  The  next  proposition  can  easily  be  proved  (see 
(20)).  It  gives  an  explicit  formula  for  the  polar  cone  of  a finitely  generated  cone. 

Pro^xisition  2.1.  Let  G £ . Then 

* n T 

pos(G)  = {x  e R ; X G ^ 0}  . 

The  following  observation  is  important.  Namely,  if  the  matrix  G e has  linearly 

* 

independent  columns,  then  the  cone  pos(G)  is  pointed  and  its  polar  cone  pos(G)  must  have 
a nonempty  interior.  This  latter  property  is  an  immediate  consequence  of  a standaid 
alternative  theorem  on  the  solvability  of  a homogeneous  system  of  linear  equations.  See 
[20]  e.g.  In  particular,  a simplicial  cone  is  pointed  and  its  polar  cone  has  a nonempty 
interior . 

Wo  now  return  to  orderings.  The  usual  ordering  < of  R*^  is  defined  "component- 
wise" and  determines  a "cone  of  nonnogati vity"  C - {x  c r":  0 < x)  which,  in  this  case, 
is  precif;ly  the  nonnegativo  orthant  r"  . This  ordering  is  generalized  in  the  following 
manner.  (See  Kransnoselski  [12].)  Let  C be  an  arbitrary  pointed  cone  in  r"  . The 
cone  ordering  (induced  by  the  pointed  cone  C ) is  defined  as  follows:  for  x,yr  r"  , 

X y if  and  only  if  y - x c C . It  can  readily  be  verified  that  this  is  a partial 
ordering.  Indeed,  axiom  (PI)  follows  from  (Cl),  (P2)  from  (C3)  and  (P3)  from  the  point- 
edness of  C . Under  this  cone  ordering  C becomes  the  cone  of  nonnegativity,  i.e. 

C = {x  c r":  0 x}  . It  is  clear  t)iat  the  cone  ordering  induced  by  the  nonnegative 

orth.int  is  precisely  the  usual  ordering.  If  C is  the  polyhedral  cone  {x  c r":  Fx  ^ O} 
where  F £ R*"^"  and  the  columns  of  F are  linearly  independent  (in  fact,  the  linear  in- 
dependence is  equivalent  to  the  pointednoss  of  C ) , then  the  cone  ordering  can  be 

viewed  as  a replica  of  tlie  usual  ordering  as  applied  to  a transformed  image  of  r"  . In- 
deed, consider  the  linear  transformation  L:  r”  -►  r""  dcined  by  L(x)  = Fx  . Then  the 
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ordering  is  precisely  the  "inverse  image"  of  the  usual  ordering  on  the  image  space 

L(r");  in  other  words,  x y if  and  only  if  L(x)  £ L(y)  . 

The  idea  of  incorporating  cone  orderings  in  the  study  of  complementarity  problems 
is  no  new  subject  at  all.  In  fact,  the  generalized  complementarity  problem  which  was  in- 
troduced by  Hcibetler  and  Price  [9]  and  later  refined  by  Karamardian  [10]  is  defined  by 
means  of  a cone  ordering.  We  shall  not  discuss  this  latter  problem  further  but  refer  the 
interested  reader  to  the  paper  [10] . 

2.2.  Known  Results.  In  this  subsection,  we  review  some  known  results  that  are  of 
fundamental  importance  in  this  paper.  In  each  of  the  four  theorems  stated  below,  the  word 
"least"  is  meant  "least  under  the  usual  ordering",  i.e.  "least  componentwise".  The  first 
theorem  characterizes  polyhedral  sets  having  least  elements. 

Theorem  2.2.  (Cottle  and  Veinott  [6])  Let  A € R*^^"  and  = {x  e r":  Ax  > b}  for 

b £ R . The  following  are  equivalent; 

(2.1a)  has  a least  element  for  each  b such  that  X^^  is  nonempty. 

T — 1 

(2.1b)  There  is  a basis  B in  A for  which  B c £ 0 for  some  c > 0 and  each  such 
basis  has  a nonnegative  inverse  (entrywise) . 

It  was  noted  in  [6]  that  X^^  is  nonempty  for  all  b if  and  only  if  there  is  an  x 
such  that  Ax  > 0 . Applying  Theorem  2.2  to  the  linear  complementarity  problem  (q,M)  , 
we  obtain  the  following  characterization  of  K-matrices. 

Theorem  2.3.  (Cottle  and  Veinott  [6])  Let  M e r"^”  . The  following  are  equivalent: 
(2.2a)  M is  a K-matrix. 

(2.2b)  For  each  q c R , X(q,M)  has  a least  element  x and  x is  the  unique  element 

— T _ 

in  X(q,M)  satisfying  x (q  + Mx)  = 0 . 

A similar  characterization  of  Z-ma trices  is  given  by 

Theorem  2.4.  (Tamir  [21])  Let  H £ R*^"  . The  following  are  equivalent: 

(2.3a)  M is  a Z-matrix. 

(2.3b)  For  each  q e r"  such  that  X(q,M)  (|i  , X(q,M)  has  a least  element  x and 

— — T — 

X satisfies  x (q  -t  Mx)  = 0 . 


.f 

I 

"V 


The  phrase  "linear  complementarity  problems  having  least-element  solutions"  used  in 
the  introduction,  originates  from  these  last  two  theorems.  It  is  an  abbreviation 
for  "linear  complementarity  problems  which  have  solutions  that  can  be  characterized  as  the 
least  elements  of  the  feasible  sets  of  the  problems." 

The  class  C of  matrices  consists  of  those  real  nXn  matrices  M which,  together 
with  Z-matrices  X and  Y,  satisfy  the  following  two  conditions; 

(Ml)  MX  = Y 

(M2)  r’'x  + s"^Y  > 0 for  some  r , s t . 

+ 

Proposition  2.5.  (Cottle  and  Pang  [4])  Let  M f C • Let  X and  Y be  Z-matrices 

satisfying  conditions  (Ml)  and  (M2) . Then 

(2.4a)  X is  nonsingular 

T T 

(2.4b)  (X  , Y ) is  a Loontief  matrix. 

Theorem  2.6.  (Cottle  and  Pang  [4])  I.et  M,  X and  Y satisfy  the  assumptions  in  Propo- 
sition 2.5.  Suppose  that  the  linear  complementarity  problem  (q,M)  is  feasible.  Then 
the  polyhedral  set 

(2.5)  '^q  “ i O'  q + Yv  2 O] 

has  a least  element  v . Furthermore  tlie  vector  x = Xv  solves  the  problem  (q,M)  and 
X can  bo  obtained  by  solving  the  linear  program 

T 

(2-6)  minimize  p x subject  to  x e X(q,M) 

T 

where  p satisfies  p X > 0 . 
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3.  POLYHEDRAL  SETS  HAVING  LEAST  ELEMENTS 


From  now  on,  we  consider  R as  being  partially  ordered  by  some  general  partial 
ordering.  We  first  state  a proposition  which  shows  how  the  least  element  of  a set  S r” 
can  be  obtained,  provided  that  it  exists.  The  proof  of  the  proposition  is  straightforward 
and  thus  omitted. 

Proposition  3.1.  Let  < be  a partial  ordering  on  R*'  . Let  x be  the  least  element  of 
a set  S £ (r",  < ) . Then  x is  a solution  to  the  minimization  problem 

(3.1)  minimize  f(x)  subject  to  x e S 

for  every  isotone  real-valued  function  f defined  on  S . Fu:thermore,  x is  the  unique 
solution  if  f is  strictly  isotone. 

The  significance  of  the  proposition  is  that  it  suggests  a con.,tructive  (though 
sometimes  not  too  effective)  approach  to  find  the  least  element  whenever  it  exists.  If  it 
happens  that  S is  a polyhedral  set  in  R*^,  and  if  f is  cliosen  to  be  linear,  tlien 
problem  (3.1)  reduces  to  a linear  program  which  of  course,  can  be  solved  by  the  simplex 
method  of  linear  programming. 

In  what  follows,  we  present  a representation  theorem  for  linear  (strictly)  isotone 
real-valued  function  defined  on  (R^,  <^)  where  C is  some  finitely  generated  pointed 


Theorem  3.2.  Let  G e have  linearly  independent  columns  and  let  C = pos(G)  . 

Consider  the  poset  (r",  ^ ^ linear  function  f:  R^'  R is  (strictly)  isotone  if 

* * »p 

and  only  if  there  exists  a vector  p e (int  C ) C such  tliat  f(x)  = p x for  every 
X € r"  . 

Proof : We  prove  only  the  case  of  strict  isotonicity.  Proposition  2.1  implies  that 

int  C = {x  £ r":  x'^'g  > 0}  . 

T * r. 

If  f(x)  = p X for  some  p £ int  C and  for  every  x £ R , then 


X <^y*»y-x  = Gr  for  some  r £ r"  \ (o) 

=»  f(y)  - f(x)  = p’’(y  - x)  - (p’^Or  > 0 . 
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Thus 


X < ^ y ■»  f (x)  < f (y)  . 

Hence  f is  strictly  isotone.  Conversely,  suppose  that  f is  strictly  isotone.  The 

linearity  of  f implies  the  existence  of  a vector  p e r"  such  that  f(x)  = p"^x  for 

every  x e r"  . It  remains  to  verify  pc  int  C . Let  r c r'J|  \ {o}  . Then  0 Gr 

because  G has  linearly  independent  columns.  The  strict  isotonicity  of  f therefore 

T 

implies  that  f(Gr)  = (p  G)r  > 0 . In  other  words,  wo  have  proved  the  implication 

r c r|J|  \ (o)  =►  (p^'or  > 0 

T * 

which  IS  clearly  equivalent  to  p G > 0 . Therefore  p c int  C . This  completes  the 
proof  of  the  theorem. 

Combining  Theorem  3.2  with  Proposition  3.1,  we  obtain 
Coroll  ai-y  3,3.  Lot  C be  given  in  Theorem  3.2.  Lot  x bo  the  least  clement  of  Sc 


(R  » ^ ^ is  a solution  to  the  minimization  problem 

...  T 

minimize  p x subject  to  x c S 

* * 

for  every  p c C . Purthermorc,  it  is  the  unique  solution  if  p c int  C 

Remark . In  both  Theorem  3.2  and  Corollary  3.3,  the  cone  C is  not  required  to  be 


simplicial. 

The  conclusion  in  Theorem  2.6  about  how  the  vector  x can  be  obtained  is  an  im- 
mediate consequence  of  Corollary  3.3  because  as  wo  shall  see  later,  x is  indct:d  the 
least  clement  of  the  feasible  set  X(q,M)  under  the  cone  order  < pQ^(X)  ' 

In  order  to  charnctori  zc  polyhedral  sets  having  least  elements  with  respect  to  cone 
orderings,  we  introduce  the  following  two  definitions. 

Definition  3.4.  (Saigal  [19])  Let  C be  a convex  set  in  r"  with  int  C y‘  <}■  . Let 
n X e 

A e R with  full  row  rank.  We  say  that  A has  the  generalized  I/ionticf  property 


with  respect  to  C if  the  following  two  conditions  are  satisfied: 
(3.2a)  there  is  a basis  B of  A such  that  C c pos(B)  j 


(3.2b)  for  each  basis  B of  A such  that  int  C n pos(B)  }<  ifi  , we  have  C £ pos(B) 


I 
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We  define  L(C)  to  be  the  set  of  all  matrices  which  have  the  generalized  Leontief  proper- 
ty with  respect  to  C . 

Definition  3.5.  (Saigal  [19])  We  say  that  a matrix  A € L(C)  is  hidden  Leontief  if 
there  is  a nonsingular  matrix  D such  that  DA  is  Leontief  and  C £ pos (D  . 

Motivation  to  study  the  generalized  Leontief  property  and  hidden  Leontief  matrices 
is  due  to  the  fact  that  "there  are  constraint  sets  {x  e Ax  = bl  that  do  not  appear 

to  be  Leontief  substitution  systems  but  can  be  shown  to  be  equivalent  to  such  systems". 
Various  characterizations  of  hidden  Leontief  matrices  have  been  obtained  by  Saigal  (191. 
See  also  [11].  An  application  which  exploits  hidden  Leontief  properties  is  given  in  [11]. 
If  A is  hidden  Leontief,  then  the  system 

Ax  = b,  X ^ 0 

is  called  a hidden  I<?ontief  substitution  system. 

Theorem  3.6.  Let  C be  a simpllcial  cone  in  . Consider  (R^,  . Let 

A f and  ~ {x  f r'':  Ax  ^ b)  for  b £ R^  . Then  the  following  are  equiv- 

alent : 

n.Ja)  has  a least  element  for  each  b such  that  ^ (J)  . 

(3. lb)  A^  € L(C*)  . 

I 

Proof : For  each  b r R , let 

{ y £ r"  : AXy  > b } 

b — 

whei'  C “ pos(X)  with  X £ R*’ ^ " and  nonsingular.  It  is  clear  that  X.  is  nonempty 

b 

<f  a id  only  if  is  so;  and  x is  the  least  element  of  X^^  with  respect  to  < ^ if 

and  only  if  y » X ^x  is  the  least  element  of  Y with  respect  to  the  usual  ordering. 

b 

Hence,  (3.3a)  is  equivalent  to 

(3.4)  for  each  b such  that  Y,  ^ 4i  , Y.  has  a least  element  with  respect  to  the  usual 

b b 

ordering. 

According  to  Theorem  2.2,  (3.4)  is  equivalent  to 

T -1 

(3.5)  there  is  a basis  B'  in  (AX)  for  which  (B')  c _>  0 for  some  c > 0 and  each 
such  basis  B'  has  a nonnegative  inverse. 
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Every  basis  B'  of  (AX)'*'  has  the  form  B'  = (BX)"*"  where  b"*"  is  a basis  of 

T T T 

The  converse  is  also  true,  i.e.,  if  B is  a basis  of  A , then  B'  = (BX)  is  a basis 
T 

of  (AX).  . Therefore,  (3.5)  is  equivalent  to 

T T T T 

(3.6)  there  is  a basis  B of  A for  which  (B  d)  X > 0 for  some  d > 0 and  each 

T -T 

such  basis  B satisfies  [c  ^ 0 =»  (BX)  c ^ 0]  . 

Noting  that 

(BX)  "*"0  ^ 0 X "*^c  = B"*^d  for  some  d e 

♦ -"P 

and  that  C = pos(X  ),  we  conclude  readily  that  (3.6)  is  indeed  equivalent  to  (3.3b), 
This  completes  the  proof  of  the  theorem. 

Theorem  3.6  characterizes  polyhedral  sets  having  least  elements  with  respect  to 
partial  orderings  induced  by  simiJlicial  cones.  It  generalizes  Theorem  2.2.  It  should 
be  fxjintcd  out  tliat  t)ie  requirement  that  tlie  cone  C be  simplicial  is  essential  in  order 

for  the  one-to-one  correspondener  hetwei  n elements  in  and  and  also  for  the 

T T 

relationship  ))etw<  cii  basis  of  (AX)  and  of  A to  lie  valid.  If  C is  merely  finitely 
goneratfd  and  pointed,  we  have  the  following  result. 

Proposition  3.7.  I,et  G c ^ ™ have  linearly  independent  columns  and  C = pos(G)  . 
Consider  (r"^,  <^)  . het  A t R^*'  . If  A**"  c L(C  ) and  is  hidden  heontief,  then 
(3.3a)  holds. 

Proof:  According  to  a property  of  hidden  heontief  matrices  {191,  there  exists  a simplicial 

r n , T * * * 

cone  S = IX  c R Dx  > 0}  such  that  A e L(.5)  and  C £ S . Since  S = (S  ) (see 

* 

[20]),  it  follows  from  Theorem  3.6  that  (3.3a)  holds  with  respect  to  S . Therefore, 

* * * 

(3.3a)  must  hold  with  re.spect  to  C because  S £ (C  ) = C . This  completes  the  proof 


of  the  proposition. 


4.  CHARACTERIZATION  OF  C 


In  this  section,  we  use  the  results  developed  in  the  last  section  to  establish  a 


characterization  of  C in  terms  of  linear  complementarity  problems  having  least-element 


solutions.  Before  proving  the  main  theorem,  we  state  and  prove  the  following  proposition 


which  describes  a relationship  between  matrices  in  C and  hidden  Leontief  matrices.  The 


proposition  generalizes  the  fact  that  if  M is  a Z-matrix,  then  the  matrix  (I,  M ) is 


Leontief. 


Proposition  4.1.  Let  M £ R n C . Then  there  exists  a simplicial  cone  C such 


T T • 

that  the  matrix  A = (I,  M ) e L(C  ) and  is  hidden  Leontief. 


Proof ; I.et  X and  Y be  Z-matrices  satisfying  conditions  (Ml)  and  (M2)  . Observe 


that  (M2)  can  be  written  as 


T T n 

(r  + s M)X  > 0 for  some  r,s  £ . 


Let  C = pos (X)  . Then  C = pos(X  ) is  simplicial.  According  to  Theorem  2.6, 


for  every  vector  q such  that  X(q,M)  ^ 4>,  the  polyhedral  set  V 


defined  by  (2.5)  has  a least  element  v with  respect  to  the  usual  ordering.  As  mentioned 


in  Theorem  3.6,  v is  the  least  element  of  V with  respect  to  the  usual  ordering  if  and 

q 


only  if  X = XV  is  the  least  element  of  X(q,M)  with  respect  to  the  cone  ordering  . 


T ’ * T 

Therefore  by  Theorem  3.6,  we  conclude  that  A £ L(C  ) . It  remains  to  verify  that  A 


IJ*  IJ*  ^ IJt 

is  hidden  Leontief.  Noting  that  (X  , Y ) = X A and  letting  D = X , we  deduce,  by 


T * -1 

(2.4b),  that  DA  is  Leontief.  Finally,  it  is  clear  that  c = pos(D  ) . This  completes 


the  proof  of  the  proposition. 


T T 

Remark;  The  matrix  A = (I,M  ) arises  in  the  linear  programming  formulation  of  the 


linear  complementarity  problem  (q,M)  with  M eC  . Indeed,  the  dual  of  the  linear  pro- 


gram (2.6)  is  given  by: 


T T 

minimize  q y subject  to  p-My^O,  y^O 


or  equivalently. 


minimize  (°)  (^)  subject  to  A*^(^)  « p,  (^)  ^ 0 


I 


Combining  Proposition  4.1  andTheorem  2.6,  wo  conclude  that  if  M eg  if  the  linear  comple- 

mentarity problem  (q,M)  is  feasible,  then  (q,M)  has  a solution  which  can  be  obtained 
by  solving  a linear  program  whose  dual  has  a constraint  set  defined  by  a hidden  Leontief 
substitution  system. 

Theorem  4.2.  Let  M e R . The  following  are  equivalent; 

(4.1a)  M e C . 

(4.1b)  There  exists  a simplicial  cone  C,  such  that  for  each  q f r"  for  which 
X(q,M)  ^ i{>,  X(q,M)  has  a least  element  x with  respect  to  < 

C 

and  X satisfies  x ”^(q  + Mx)  = 0 . 


Proof:  (4.1a)  =*  (4.1b).  This  follows  immediately  from  Theorem  2.6  and  the  proof  of 

Proposition  4.1. 

(4.1b)  (4.1a)  . Lot  C = pos(X)  where  X £ is  nonsingular  but  not  neces- 


sarily  a Z- 

•matrix.  Let 

Y - MX  . According  to  the  assumption,  we  deduce  tliat  for  every 

vector  q 

for  which  the 

set  V (defined  by  (2.5)) 

q 

is  nonempty, 

, it  has  a least  element 

-1_ 

v(=  X X) 

with  respect 

to  the  usual  ordering  and  v 

satisfies 

(Xv)’^  (q  + Yv)  = 0 . 

I<et  )c  bo 

an  index  in 

{l,...,n}  . Choose  q = e*^  - 

- Ya’' 

whore 

k -Ik 

a = X e . Then 

)c 

a £ V . 
q 

Tliereforc  V 

q 

. — k 

contains  an  element  v < 

k 

a 

k k 

Moreover,  v a . Define 

k k - 

-k  k 

V = a - \ 

r , then  v 

is  a non-vanis)iing  vector. 

For 

•H 

, we  have 

(4.2a) 

X.  v'^  = X.  (x'^e*^)  - X.  v’' 

< 0 

!•  1*  i« 

and 

(4.2b) 

Y.  v^  = Y.  a*^  - Y.  v’' 

= q.  - Y v'^  < 0 . 

1 1*  — 


• • Ink 

Now  if  wo  define  the  matrix  w = (v  ) where  v are  the  vectors  defined  above, 

then  clearly,  X'  = XW  and  Y*  = YW  are  Z-matrices  by  (4.2),  Moreover  MX*  = Y*  . It 

remains  to  verify  that  there  exist  vectors  r and  s in  r"  such  that  r^X*  + s^Y'  > 0 . 

By  Theorem  3.6,  it  follows  that  = (I,  m"^)  e l(C*)  . Therefore,  there  exists  a basis 

T T * T * 

B of  A such  that  C £ pos(B  ) . This  implies,  by  the  fact  that  int  C ^ 4*  . 

(t'"  B)X  > 0 for  some  t f rJJ[ 
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If  we  define  the  vectors  r = (r^)  and  s = (s^) 


* i T 

t.  if  e *=  (B  ) . for  some  j 

3 *3 


otherwise 


t.  if  (m"^)  . = (b"^)  for  some  j 

3 •!  *3 


0 otherwise. 


then  clearly  r,  s e and 


Finally,  we  have 


w ip  7 *r 

rX  + sY=(r  + s M)X  = (t  B)X  > 0 


T T T T 

r X’  + s y = (r  X + s Y)W  > 0 


because  each  column  of  W is  non-vanishing  and  W ^ 0 . This  shows  that  M e C and 
completes  the  proof  of  the  theorem. 

The  above  theorem  generalizes  both  Theorems  2.3  and  2.4.  It  should  be  emphasized 
that  in  (4.1b),  the  cone  C is  not  required  to  be  induced  by  a Z-matrix.  Moreover,  it  is 
worth  pointing  out  that  there  exist  matrices  M e C which  satisfy  the  defining  conditions 
(Ml)  and  (m2)  for  some  matrices  X and  Y which  are  not  both  Z-matrices.  An  example  is 
given  by  the  following. 

Example.  Consider  M = (^  ^)  . Let  X = ( ,^  “h  and  y = ( ? ~\)  ■ Clearly 

u -1  -/j  1 /j  -1 

MX  = Y and  (M2)  is  satisfied  because  X is  a K-matrix.  M e C because 


^)  (^  ~^)  = (^  ”^) 

'o  -l'  'o  l'  'o  -l' 


Observe  that  Y is  not  a Z-matrix. 
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1 20.  ABSTRACT  fConrinuo  on  ofdo  If  nacoaosry  «nd  Idenilfy  by  6iorA  mimbor; 

I ^This  paper  first  generalizes  a characterization  of  polyhedral  sets  having 
least  elements,  wliich  is  obtained  by  Cottle  and  Veinott  to  the  situation 

wlicre  Euclidean  space  is  partially  ordered  by  some  general  cone  ordering  (rather 
than  the  usual  ordering).  We  then  use  this  generalization  to  csteiblish  the  fol- 
lowing characterization  of  the  class  C of  matrices  ( C arises, as  a generaliza- 
tion of  the  class  of  Z-matriccs,  •ce  Hli  l-Hlr-L1414t  1^1  c C'y^f  and  only  if  foij 
every  vector  q for  which  the  linear  complementarity  problem  (q,M)  is  feasible 
the  problem  (q.H)  has  a solution  which  is  tho  least  element  of  the  feasible  set 
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of  (q,M)  with  respect  to  a cone  ordering  induced  by  some  simplicial  cone.  This 
latter  result  generalizes  the  characterizations  of  K-  and  2-matriccs  obtained  by 
Cottle  and  Veinott  [6J  and  Tamir  [21]  respectively. 
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